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ABSTRACT 

Microlensing light curves are now being monitored with the temporal sampling and photometric preci- 
sion required to detect small perturbations due to planetary companions of the primary lens. Microlens- 
ing is complementary to other planetary search techniques, both in the mass and orbital separation of the 
planets to which it is sensitive and its potential for measuring the statistical frequency of planets beyond 
the solar neighborhood. We present an algorithm to analyze the efficiency with which the presence of 
lensing binaries of given mass ratio and angular separation can be detected in real microlensing datasets. 
Such an analysis is required in order to draw statistical inferences about lensing companions, and differs 
from previous studies of idealized microlensing experiments by incorporating instead the actual sam- 
pling, photometric precision, and monitored duration of individual light curves. We apply the method 
to artificial (but realistic) data to explore the dependence of detection efficiencies on observational pa- 
rameters, the impact parameter of the event, the finite size of the background source, the amount of 
unlensed (blended) light, and the criterion used to define a detection. We find that: (1) the integrated 
efficiency depends strongly on the impact parameter distribution of the monitored events, (2) calculated 
detection efficiencies are robust to changes in detection criterion for strict criteria (Ax^ ~ 100) and large 
mass ratios {q~ 10~^), (3) finite sources can dramatically alter detection efficiencies to companions with 
mass ratios q~ 10~^, and (4) accurate determination of the blended light fraction is crucial for the accu- 
rate determination of the detection efficiency of individual events. Suggestions are given for addressing 
complications associated with computing accurate detection efficiencies of real datasets. 

Subject headings: planetary systems — binaries — gravitational lensing 



1. INTRODUCTION 

Th e discovery in 1995 of a massive planet orbiting 51 
' i3i 



P e g ( M ayor fc Qucloz 19 95 ) followed by the discovery of 



oovoral more planoto orbiting noarb y dwarf otaro uoing the 



have been made, but the measurements are difficult and 
the detections remain controversial; planned s pace-based 
missions are expe cted to be more successful (Lindegren 
fc Ferryman 1996). Occultation methods use very accu- 



same radial velocity technique (see [Marcy fc Butler 1998 
and references therein) has focussed both public and sci- 
entific attention on the search for extra-solar planets and 
the experimental and theoretical progress being made in 
developing other viable detection techniques. 

Duo to their omall maoo and DiHO, extra oolar planoto are 



rate photometry of the parent star to detect the small 
decrease in flux (5 1%) caused by a planet transiting the 



difficult to find , Proposed detection methods can be subdi- 



face of the star ( Hale fc Doyle 1994). Occultation se arches 
are currently being carried out ( Deeg et al. 1998 ), and 
spaced-based missions are pl anned to increase the sens! 
tivitv to sm all-mass planets ( Deleuil et al. 1997 , Borucki 
et al. 1997). By far the most successful indirect method 



vided into direct and indirect techniques. Direct methods 
rely on the detection of the reflected light of the parent 
star, and are exceedingly challenging due to the extremely 
small flux expected from the planet, which is overwhelme d 
by stray light from the star itself ( Angel fc Woolf 1997 ). 



i Somc direct imaging searches have alrcad^^" been performed 



for discovering planets has been the Doppler technique, 
which employs extraordinarily precise radial velocity mea- 
surements of nearby stars to detect Doppler shifts caused 
by orbiting planets. Several teams have monitored nearby 
stars with the aim of detecting the D o ppler signal of or 
biting planets ( McMillan et al. 1993 \_ Mayor fc Queloz 
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199f 



the cons truction and launchi ng of a space-based interfer- 
ometer ( |Woolf fc Angel 1998| ). 

Astrometric, radial velocity and occultation measure- 
ments can be used to detect the presence of a planet in- 
directly. Astrometric detection relies on the measurement 
of the wobble of the stellar centroid caused by the motion 
of the star around the center-of-mass of the planet-star 
system, and yields the mass ratio of the planet-star sys- 
tem. Many attempts to find extra-solar planets in this way 
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Butler et al. 1996, Cochran et al. 1997, Noyes et al 
ib date these groups combined have discoverec 
over 20 extra-solar planets using this technique, with new 
planetarv companions b eing announced every few months 
( Marcy fc Butler 1998| ). Of these techniques, only the 
proposed space-based interferometric imaging and transit 
searches are expected to be sensitive to Earth-mass plan- 
ets. 

These detection techniques are complementary to one 
another both in terms of their sensitivity to planetary mass 
and orbital separations and the specific physical quantities 
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of the planetary system that they measure. All share two 
distinct advantages: the experiments are repeatable and, 
due to their reliance on flux measurements of the parent 
star or the planet itself, they are sensitive to stars in the so- 
lar neighborhood where follow-up studies can be most eas- 
ily pursued. For example, spectroscopic follow-up studies 
may enable the detection of molecules commonly thought 
to be indicative of life, such as water, carbon dioxide, and 



ozone ( Woolf fc Angel 1998 ). This advantage is linked to 
a common drawback: the searches can be conducted only 
on a limited number of nearby stars. Without substan- 
tial developments in technology beyond what is currently 
available, these methods are unlikely to detect planets or- 
biting stars more distant than ^1 kpc, and are thus unable 
to address questions about the nature of planetary systems 
beyond the immediate solar neighborhood. 

Microlensing was proposed in 1986 by Paczyhski as a 
method to detect compact baryonic dark matter in the 
halo of our Galaxy. Soon a fter, thre e collaborations (MA 



CHO, Alcock et al. 1993; EROS, Aubourg et al. 1993 



OGL E, [Udalski et al. 1993) were organized to implement 



this proposal with massive observational programs. The 



basic idea is simple: when a compact object, such as a 
star or massive dark object, passes near the line of sight 
to a distant source, the gravitational field of the interven- 
ing object will serve as a lens, creating two images of the 
distant source. For typical stellar sources and lenses in our 
Galaxy, the separation of the images will be on the order of 
milliarcseconds and hence unresolved. However, the lens 
will a lso magnify the source; this magnification is mcasur- 



able and depends (only) on the angular separation of the 
lens and source. Since the lens, source, and observer are all 
in relative motion, the magnification will be time- variable, 
creating a "microlensing event." Because the probability 
that any single source star will be microlensed (the 'op- 
tical depth') is so low, 0(10"^), millions of stars must 
be monitored; the crowded fields towards the Large and 
Small Magellenic clouds (LMC and SMC) and the Galac- 
tic bulge were thus the natural targets. In 1993, the first 
candid ate microlensing events were announced toward the 
LMC ( lAlcock ct al. 1993t [Aubourg et al. 1993D and Galac- 
tic bulge (Udalski ct al. 1993). Six years later, these three 
'survey teams' continue to search and discover microlens- 
ing events. Over 400 events have been discovered toward 
all three targ ets (MACHO, Alcock et al. 1997a, 1997b, 
1997c : ERO S, [Renault et al. 199^, Falanque-Delabrouille 



et al. |l99^ ; OGLE, [Udalski et al. 1997] ), the overwhelm- 



ing majority of which are discovered toward the Galactic 
bulge. 

Although Paczyhski's original suggestion was to search 
for dark matter, before the first event was discovered Mao 
& Paczyhski (1991) had already noted that it might be 
possible to detect planetary companions of the primary 
microlenses via the distortions they create in the magnifi- 
cation pattern, and thus the light curve, generated by the 
primary lens. The nature of the distortion depends on the 
mass ratio and angular separation of the two components. 

As a planet search technique, microlensing offers unique 
advantages. Since microlensing is caused the gravitational 
field (i.e., mass) of the lenses, it is not hmited to the study 
of nearby or indeed luminous objects, and thus can be used 
to search for planetary companions around typical Galac- 
tic stars at distances of many kiloparsecs. As a conse- 



quence, a nearly unlimited number of dwarf stars are avail- 
able to serve as gravitational microlenses and potential 
search candidates. This advantage is linked to the primary 
drawback of microlensing planetary searches: most follow- 
up studies will be difficult due to the faintness of the stars 
serving as typical lenses. This drawback is compounded 
by the irrepeatability of specific microlensing observations; 
lensing of a particular source by a particular lens is a sin- 
gular occurrence. Nevertheless, the robust statistics on the 
nature of planetary systems many kiloparsecs distant and 
the complementary nature of the information about dis- 
covered systems microlensing can provide make it an im- 
portant tool i n the cadre bei ng assembled to study extra- 
solar planets (Sackctt 1999a). Furthermore, microlensing 
planet searches are relatively inexpensive, requiring only 
several dedicated 1-meter class telescopes. Microlensing is 
the only technique currently capable of routinely discov- 
ering planets like our own Jupiter, and the only ground- 
based method capable, in principle, of detecting distant 
terrestrial-mass planets, though this will requi re substan- 
tial enhancements over existing cap abilities (Bennett & 
Rhic 1996[ ; [Pcale 1997| , [Sackett 1997| ) 



Current microlensing survey teams have sampling rates 
that are too large ('^day) and/or photometric accura- 
cies that are too poor 5%) to detect and characterize 
these perturbations, but the nearly 100 real-time electronic 
alerts of on-going bulge events that they provide annually 
have become the primary targets for newly-formed mi- 
crolensing planet searches. These new 'monit oring teams 



(PLAN ET, A lbrow et al. 1997 1998; OMAN, [Alcock et al 
1997d[ ; MPS, [Rhie et al. 19991 ) have formed with the ex- 



press purpose of executing the nearly continuous temporal 
coverage and high photometric precision on real-time mi- 
crolensing alerts necessary to detect deviations from the 
generic light curve of the sort expecte d from planetary an d 
other microlensing anomalies (see e.g. Albrow et al. 1998). 
The PLANET collaboration in particular has now moni- 
tored nearly 100 microlensing events with varying degrees 
of photometric sampling and precision (often '^hourly with 
^2% over the largest magnification regions); such existing 
datasets may already place interesting constraints on the 
frequency of stellar binaries and planetary systems. 

Like all planet search techniques, microlensing is not 
100% efficient, due to both intrinsic and observational lim- 
itations. The efficiency with which a given dataset will re- 
veal the presence of a companion to the primary microlens 
must be quantified before it can be used to constrain the 
frequency and properties of extra-solar planetary systems. 
Quantification of detection efficiencies of any kind can be 
a difficult and tedious process: the intrinsic limitations 
of the method must be identified and combined with the 
actual observational limitations. The detection efficiency 
may depend on hidden or unmeasurable parameters; these 
must be identified and properly quantified in order to avoid 
biasing the final conclusions. Detection efficiencies of other 
planet se arch techniques have b een presented and applied 
to data (Nelson & Angel 1998). However, despite a sub- 
stantial body of work addressing the likely planet detec- 
tion efficiency of idealized microlensing programs (see § 
3), no methods have yet been proposed for calculating the 
efficiency of actual datasets to lensing binaries of a given 
type. 

Here we present an algorithm for computing the detec- 
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tion efRciency to lensing binaries that is specific to indi- 
vidual microlensing light curves. By directly imposing the 
actual observational limitations as a constraint, the ap- 
proach is less prone to the biases that may arise when 
using simplified models of observational conditions. Fur- 
thermore, the method is simple to implement and compu- 
tationally inexpensive, since it involves direct integration 
over unknown quantities rather than Monte Carlo simula- 
tions commonly used to calculate detection efficiencies of 
idealized observing programs. We apply this method to 
simulated datasets in order to explore how the detection 
efRciency depends on intrinsic and observational effects. 
We also explore possible biases that may be introduced 
into the inferred efficiency of individual events to planet 
detection if the size of the source and the fraction of un- 
resolved light ('blending') are ill-constrained. 

In § 2, we review the relevant formalism for microlens- 
ing by single and double stars. A brief review of the lit- 
erature on expectations for idealized microlensing planet 
searches is presented in § 3, contrasting these to the goal 
of this paper. Detection and detection efficiency, as used 
throughout this work, are defined in § 4, along with a 
general overview of the connection between detection effi- 
ciency and the mass ratio, angular separation and impact 
parameter of the event. Our algorithm for computing de- 
tection efficiencies is described in § 5 and applied to artifi- 
cial data in order to access the effects of different detection 
criterion. The effects of finite source size on detection ef- 
ficiencies are presented in § 6; the effects of blending in 
§ 7. Suggestions for addressing complications associated 
with computing accurate detection efficiencies of real mi- 
crolensing datasets are given in § 8. We summarize and 
conclude in § 9. 

2. RELEVANT FORMALISM FOR SINGLE AND DOUBLE 
MICROLENSES 

2.1. Single lenses 

The time- variable flux observed from a microlensed star 



IS 



F{t) ^ Fo[Ait) + fE 



(1) 



where Fq is the unlensed flux of the star, /s is the ratio 
of any unresolved, unlensed background light to Fq (the 
'blend fraction'), and A{t) is the magnification. The mag- 
nification of Aq of a point source by a point lens can be 
written as 

u{u^ + 4)1/2 

Here u is the instantaneous angular separation of the 
source and the lens in units of the angular Einstein ring 
radius Oe oi the lens, a degenerate combination of the lens 
mass and and distances defined by 



4GM Dls 



F>oL F>os 



1/2 



130 /ias 



M 



1/2 



0.2Mq 

(3) 

where M is the mass of the lens, and I?ls, F>os, Dql 
are the lens-source, observer-source, and observer-lens dis- 
tances, respectively. For the scaling relation on the far 
right of equation (g) , we have assumed Dql = 6 kpc and 
Dos = 8 kpc. Note that for u ^ 0, Aq 1/u. Since the 
source, lens, and observer are all in relative motion, u will 



be a function of time. For rectilinear motion, it can be 
shown that 



u{t) = 



tE 



-I 1/2 



(4) 



where to is the time of maximum magnification, Umin is 
the minimum angular separation, or impact parameter, of 
the event in units oi 9e, and Ie is a characteristic time 
scale of the event, the Einstein time, defined by 



SeDql 



10 days 



M 



0.2M, 



1/2 



(5) 



Here v± is the transverse velocity of the lens relative to 
the observer-source line-of-sight. For the scaling rela- 
tion on the far right of equation (|^), we have assumed 
Doh = 6 kpc, Dos = 8 kpc, and v± = 100 kms"^. 

A point-lcns point-source (PSPL) light curve is thus a 
function of five parameters, to, tE, Umin, ^o, and fs- Un- 
less the lens is measurably luminous, the parameters Fq 
and /b depend only on the source and its environment. 
The distributions of Fq and Jb depend on the luminosity 
functions of the observed sources and any blended light. 
The parameters and Umin are purely geometrical. Since 
the distribution of to is flat and its value has no effect on 
the analysis of detection efficiencies, we will hereafter set 
to — 0. The intrinsic distribution of Wmin is also flat, but 
the observed distribution has an upper limit set by the 
detection threshold of the survey teams. If the lens is not 
contributing significantly to the blended light, only the 
characteristic time tE contains physical information about 
the lens itself. Its intrinsic distribution is set by the lens 
masses and spatial distribution of the lenses and sources; 
the observed distribution oi tE depends on the temporal 
sampling of the microlensing campaigns. An example of 
a point-lens light curve with to = 15 days, tE = 30 days, 
Umin = 0.5, Fq = 1, and /b — is shown as the dotted 
line in the top panel of Figure |l|. 

2.2. Binary lenses 

A companion orbiting a lens constitutes a double lens 
and thus can be described by the formalism of binary 
lenses. The flux is still expressed by equation (^, but 
the magnification can no longer be calculated analytically. 
Instead, the lens equation describing the mapping from 
the source plane (^,?7) to the image plane {x,y) must be 
solved. Following Witt (1990), we write the lens equation 
for two masses with fractional mass mi and m2 located at 
positions zi and Z2 in terms of the complex coordinates 
( = 7] + i^ and z = x + iy, 



mi 



zi 



+ 



m2 



Z2 



(6) 



Here all distances are in units of 6e, the Einstein ring 
radius of the total mass M of the binary. The mapping 
is described completely by two parameters, the instanta- 
neous angular separation of the two components in units 
oi 6e, b = \z\ — z-2\, and the mass ratio of the system, 
1 — m2/rni. Without loss of generality, we will assume 
that m2 < mi, so that q < 1. Equation (^ is equiva- 
lent to a fifth-order complex polynomial in z, which can 
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Fig. 1. — Top panel: Magnification as a function of time in days for a microlonsing event with (solid curve) and without (dotted curve) the 
presence of a companion with mass ratio q = 10~^. The event has an Einstein time scale of = 30 days and minimum impact parameter in 
units of the Einstein ring radius of ttmin = 0.5. The single lens curve reaches a maximum magnification of Amax = 2.2 at a time to = 15 days. 
The binary lens curve was generated assuming that the source trajectory makes angle = 300° with respect to the binary axis. Also shown 
is the best single-lens fit to the binary- lens light curve (dashed curve). Bottom panel: The fractional deviation 5 of the binary-lens light curve 
from the single-lens light curve with parameters = 30 days,Mniin = 0.5, to = 15 days (solid curve) and from the best fit single-lens light 
curve (dashed curve). 



be solved by the usual techniques. Each source position 
produces either three or five images. The magnification 
Aj of each image j is inversely proportional to the deter- 
minant of the Jacobian of the lens mapping, evaluated at 
that image position. 



1 



~ \detJ\ 



detJ.l-p|. (7) 

oz oz 



The total magnification is given by the sum of the in- 
dividual magnifications, A = J^j^j- The set of source 
positions for which the magnification is formally infinite, 
given by the condition det J = 0, defines a set of closed 
curves called caustics. Depending on the values of b and 
q, a binary lens may have one, two, or three caustics; the 
multiplicity of images changes by two as the source crosses 
a caustic. 

A static, point-source binary-lens light curve is a func- 
tion of eight parameters. Two are identical to the point- 
lens case, Fq and fs- The equation for retains the same 
form, but the choice of the fiducial mass M is arbitrary 
and can refer to the total mass of the binary or the mass 
of one of the components. For binaries, Wmin refers to the 
minimum angular separation (in units of 9e) between the 
source and the origin of the binary system. The choice of 
origin is also arbitrary; popular choices are the position 
of center of mass, the position of one of the masses, or 
the midpoint between the two. Equation (Q) still holds for 



static binaries, so that to is the time at which u = Umin, 
but for binary lenses this need not be the time of max- 
imum magnification. The mapping parameters b and g, 
and the angle 9 on the sky between source trajectory and 
the binary axis, are the final three binary lens parameters. 
In the case of a single lens, the lensing geometry is az- 
imuthally symmetric, and 6 is completely degenerate for 
any measured light curve. 

The value of Umin has a large effect on the detection effi- 
ciency of a given light curve to lens binarity; smaller Mmin 
events generally have higher efficiency. The Einstein time 
scale ts affects the detection efficiencies in that shorter 
time scale events will, in general, be less densely sampled 
by monitoring teams than longer time scale events. Blend- 
ing complicates matters due to the ambiguity between 
light curves with different combinations of itmin and blend 
fraction Jb (Wozniak & Pacyhski 1997), which can lead to 
ill-determined detection efficiencies if the blend fraction is 
poorly constrained (§ 7). 

The magnification patterns of close binaries (6 ^ 1), 
wide binaries (6^1), and binaries with small mass ratios 
{q <^ 1), can be written as mathematica l perturbations to 
the single lens pattern ( Dominik 1999 ); the light curves 
they produce can be mistaken for those due to a single 
lens for a majority of source trajectories. The fractional 
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deviation 6, which is defined by 

^^^^ ^ Ait) ^ Mt) 



A{t) 



(8) 



where A(t) is the binary-lens magnification and Ao{t) is 
the magnification of the best-fit single-lens model, quanti- 
fies the degree to which a binary-lens light curve deviates 
from a best-fit single lens model as a function of time. 

The best-fit single lens model need not have the same 
parameters Iq, Ie and Mmin as the underlying binary and 
will vary depending on the size and duration of the devi- 
ation, which is determined by the source trajectory 6 and 
the binary-lens parameters (6, q) . The top panel of Fig- 
ure shows a light curve for binary system with q — 0.01 
and b = 1.0, and two single-lens light curves: one assum- 
ing that the parameters tE,to, Umin are the same as in the 
underlying system, and the other the best-fit single lens 
curve. In this particular instance, the diff'erence between 
the two single-lens curves is not large, although this is not 
universally true (see § 5.1). The bottom panel shows the 
fractional deviation of the binary light curve from both 
single-lens light curves. In this example, the deviation is 
appreciable {S > 1%) for a large fraction of the light curve, 
but is large {6 > 5%) for only ^ 5 days, underscoring the 
need for high photometric precision in microlensing planet 
searches. 

3. DETECTION PROBABILITIES FOR IDEALIZED 
MICROLENSING PLANET SEARCHES 

The success of microlensing survey and monitoring 
teams in implementing massive observational programs 
has prompted, and in part been driven by, theoretical 
work outlining optimal observational strategies, expected 
detection probabilities for binary lens perturbations, and 
methods for characterizing and extracting additional in- 
formation about the companions. 

Using heuristic arguments based on rough scaling rela- 
tions for the size of planetary caustics, Mao & Paczynski 
(1991) estimated in their seminal paper that 3% of mi- 
crolensing light curves should cross caustics and thus show 
planetary deviations if all lens have a Jupiter-mass planet 
at instantaneous angular separations comparable to the 
Einstein ring radius of the primary. Randomizing the or- 
bital phase and inclination of the orbital plane to obtain 
the distribution of angular separations hOE on the sky cor- 
responding to a real orbital radius a, Gould & Loeb (1992) 
estimated that planets with characteristics like our own 
Jupiter and Saturn orbiting solar-type stars halfway to 
the Galactic center would have detection probabilities of 
17% and 3%, respectively. Their substantial probability 
is based on the assumption non-caustic crossing perturba- 
tions as small as 5~ 5% can be detected regardless of their 
duration. Performing a similar study, Bolatto & Falco 
(1993) obtained a higher detection probability of '^40% 
for Jovian-like planets, requiring that the integrated dif- 
ference between the binary light curve and a single-lens 
light curve of the same impact parameter Mmin exceed a 
threshold corresponding to a 10% deviation over 36 hours. 
Bennett & Rhie (1996) included for the first time the fi- 
nite size of the source, enabling them to study planetary 
systems with mass ratios as low as 10~^, comparable to 
the mass ratio of the Earth and a parent star with mass 



O.5M0, whose caustics are smaller than the angular size of 
typical source stars. They found a small but non-negligible 
detection probability of ~ 2% for such low-mass planets, 
assuming a 50% weather duty cycle for observations capa- 
ble of detecting deviations as small as 4% and as brief as 
iB/200 (~2.4 hour). 

All of these studies have concentrated on companions 
with instantaneous angular separations in the "Lousing 
Zone," a region defined by 0.6 < b < 1.6 in which plan- 
etary caustics may cross the source trajectory for events 
with Mmin < 1, corresponding to an initial alerting mag- 
nification by the survey teams of A > 1.34. Recently, Di 
Stefano & Scalzo (1999a, 1999b) have presented exhaus- 
tive studies of the probability of detecting planets in wide 
orbits that place gas giants like those in our own Solar 
System outside the Lensing Zone of typical lenses. They 
conclude that such events may be detectable as isolated 
short duration events, which will alter the time scale dis- 
tribution of all microlensing events, or as repeating events 
as the source trajectory passes first near the primary and 
then near the companion on the sky. The observing and 
analysis strategies of the survey and monitoring programs 
must be altered to optimally detect such planetary de- 
viations, which, by definition, occur when the source is 
outside the Einstein ring of the primary. 

Wambsganss (1997) has demonstrated the extraordi- 
nary variety of perturbations that can be produced by 
given planetary systems depending on the trajectory of the 
source star. In principle, this variety might lead to ambi- 
guity in the inferred properties of a lensing planet based on 
a measured light curve, a potentially serious difficulty for 
the characterization of detected planetary systems. Gaudi 
& Gould (1997) found that potentially severe fitting am- 
biguities of a factor of ^20 in the measured mass ratio 
could plague light curves with planetary anomalies, but 
that these could be alleviated by dense and accurate sam- 
pling of the planetary perturbations or simultaneous opti- 
cal/infrared photometry. Gaudi (1998) notes that plane- 
tary microlensing anomalies may also be mistaken for par- 
ticular types of binary source events; the ambiguity can 
be mitigated with similar additional observations. Griest 
& Safizadeh (1998) showed that the planetary detection 
probabilities for high magnification events can be quite 
large since the source necessarily passes close to the cen- 
tral caustic generated by any binary. Using several detec- 
tion criteria, they found that Jupiter-mass planets can be 
detected with ~100% efficiency in events with Mmin-O.l 
(^max ~ 10) over a substantial range of parent-planet sepa- 
rations. Gaudi, Naber, & Sackett (1998) pointed out that 
this necessarily implies that the presence of multiple plan- 
ets would be revealed in high magnification events, since 
all companions contribute to the central caustic. They also 
noted that unless the caustic morphology of higher multi- 
plicity lensing systems is understood, fitting degeneracies 
will prevent the association of deviant light curves with 
well-characterized planets in such events. 

Such considerations have led to detailed proposals for 
observational strategies designed to maximize the de- 
tection of planets via microlensing. For planets of all 
mass ratios, Peale (1997) has stressed the importance 
of a longitudinally-distributed network of southern tele- 
scopes for continuous monitoring, as has implemented 
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by the PLA NET collaboration si nce 1995 in its Jovian- 
mass search ( [Albrow et al. 1997| , 1998). Sackctt (1997) 
points out that the high burden of proof demanded of non- 
repeating microlensing events necessitates a high-sampling 
rate for proper characterization of planetary anomalies; 
this suggests a wide-field imager at an excellent site as an 
ideal component in any microlensing search effort aimed 
at small mass planets. Both studies conclude that tens of 
planets per year may be detected by aggressive programs 
if every microlens has a planet in its Lensing Zone. 

4. DEFINING DETECTION CRITERIA AND EFFICIENCIES 
FOR ACTUAL DATASETS 

Previous studies of detection probabilities of idealized 
microlensing planet searches (§ 3) differ in spirit and 
approach from the situation confronted by a researcher 
wishing to use an actual database of microlensing light 
curves to draw inferences about possible lensing compan- 
ions. First, real light curves are irregularly sampled, each 
with its own and usually varying photometric precision. 
Second, only the post-alert portion of the light curve is 
available to current monitoring teams, resulting in differ- 
ences in the monitored phase of each event. Third, the 
intrinsic impact parameter distribution of the events is al- 
tered by choices made by both the discovery and moni- 
toring teams; the actual distribution cannot be assumed 
to be uniform over any interval. Finally, because the true 
lensing parameters are unknown, a binary lens fit must be 
compared to the best-fit single lens model for a given light 
curve, not to a single-lens model in which the companion 
has been removed (see Fig. |^). Previous studies that did 
not use the best-fit single-lens model as the null hypoth- 
esis when calculating detection probabilities for idealized 
searches have overestimated the true detection probability 
(§ 5.1). These differences motivate the detection criteria 
and method of calculating detection efficiencies for actual 
datasets that we now describe. 

All light curves contain information about the presence 
of companions around Galactic lenses: obviously anoma- 
lous light curves signal the possible presence of a compan- 
ion while light curves without a detectable anomaly sig- 
nal the possible absence of (certain types of) companions. 
With real datasets, both statements are probabilistic: the 
presence or absence a particular type of lensing compan- 
ion in a particular system can only be made with a cer- 
tain degree of confidence. Ultimately, we seek a method 
to characterize this probability consistently for all light 
curves, whether obviously anomalous or not, so that the 
complete dataset can be used to constrain the distribution 
of planets and other companions in orbit around Galactic 
lenses. 

Because most microlensing light curves do not show ob- 
vious anomalies indicative of companions, we begin here 
by quantifying the extent to which an apparently non- 
anomalous light curve can be taken as a sign that plan- 
etary companions of a given mass ratio q and projected 
separation b are truly absent in the microlensing system 
under consideration. The absence of an observed planetary 
anomaly may be due to insufficient sampling, photometric 
precision, or monitoring duration, rather than the absence 
of a planet. This means that the efficiency of detection 
must be computed separately for each event. Intrinsic pa- 
rameters, such as source trajectory, source size, and blend- 



ing, may also serve to 'hide' the planetary anomaly from 
the observer, and must be disentangled from observational 
effects both to compute detection efficiencies of individ- 
ual datasets and to formulate future observing strategies 
which will maximize the detection probability. 

Before we can begin to quantify the detection efficiency, 
ei{b,q), of a lensing binary with microlensing parameters 
6, q using light curve £i, we must define the meaning of 'de- 
tection.' Here, we will consider a planetary companion to a 
lens to be 'detected' in light curve £i if some combination of 
the intrinsic binary parameters b, q and 9 produces a sub- 
stantially better fit (characterized by A^^) to the observed 
light curve than the best-fit single lens model. In making 
the comparison, the parameters Iq, Ie, Fq, Umin and Jb are 
allowed to vary to achieve the best fit in both the binary 
and single lens models. The meaning of 'substantially bet- 
ter fit' can be adjusted by altering the threshold value of 

^Xthrcsh that = Xbest single ~ Xbest binary mUSt exCecd. 

Since here we consider only light curves that are not ob- 
viously discrepant from PSPL (by the criterion above) , the 
angle 9 between the source trajectory and binary axis is de- 
generate, and can be assumed to be drawn from a random 
distribution over the full range Q < 9 < 2tt. By detection 
efficiency, ei{b,q), we will mean the probability that an 
actual planetary companion with mass ratio q and instan- 
taneous angular separation b would be 'detected' by the 
criterion above in a given observed light curve £i, assum- 
ing a random source trajectory. Only a subset of source 
trajectory angles will produce detections by this criterion. 
The efficiency incorporates both the intrinsic sensitivity of 
the observations to a given anomaly and the probability 
that such an anomaly occurs for random source trajecto- 
ries. An efficiency ei(b,q) of zero implies that a lensing 
companion with characteristics b and q would always es- 
cape detection with these observations; ti{b, q) ~ \ implies 
that the companions would always be detected (if present) 
with data of this type and quality regardless of source tra- 
jectory through the magnification pattern. This efficiency 
can then be used to place a confidence level on the non- 
detection of 6, q binary anomalies in a given light curve. 

Note that in order to draw inferences about the non- 
detection of planetary companions of absolute mass rup (in 
Mq) and orbital radius a (in AU), additional assumptions 
or information must be brought to bear on the distribu- 
tion of lens masses and distances. Integrations over orbital 
phase and inclination then must be performed to deproject 
the instantaneous angular separation b into a probabilistic 
distribution function for orbital radius a. 

4.1. Qualitative effect of the Intrinsic Binary Parameters 

In addition to its dependence on observational parame- 
ters, the detection efficiency of a given light curve will be 
a strong function of the intrinsic mass ratio q and instan- 
taneous binary separation 6 of the lens system. 

Three limits are of special interest. First, as 6 0, a bi- 
nary lens reduces to a single lens, so that most light curves 
will have small detection efficiencies to binaries with small 
separations. Second, for 6 ^ 1, a binary lens reduces to 
two isolated point lenses. Only small Wmin trajectories with 
~ (i.e. parallel t o the binary axis) will be tray the pres- 
ence of both lenses ( Pi Stcfano fc Mao 1996| ); most trajec- 
tories will produce light curves that are indistinguishable 
from single lens light curves. Finally, for g — > 0, a binary 
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Fig. 2. — Contours of fractional deviation <5, for binary separations b = 0.2, 0.6, 1.0, 1.5, 2.0, 3.0 and mass ratios q = 1.0, 10^^, 10^^, 10"'^, 
10~*, as a function of position in the source plane (»?, ^) in units of the Einstein ring radius of the binary, 9e- The red contours are 
5 = +0.01, +0.05 (thin and heavy lines), the blue contours are 5 = —0.01,-0.05 (thin and heavy lines), and the heavy black lines are the 
caustics, S = CO. The green circle is the Einstein ring. Black dots mark the positions of the lensing masses. The thin straight line is a 
trajectory with u^nin = 0-5 3,nd 9 = 300. Light curves for these trajectories are shown in Figure 3. 



lens reduces to a single lens. Detection efficiencies thus 
will be fiighest for g ~ 1 and b \, and will decline for 
decreasing q and outside an annulus centered near & ~ 1. 

To illustrate these points, we calculate (5 as a function 
source position , 77) for a grid of b and q values. Since the 
best-fit single lens light curve, and therefore (5, depends on 
the exact trajectory, we must adopt heuristic approxima- 
tions in order to display these qualitative results in a single 
plot of 5 contours. For & < 1.0 and all g, we choose Wmin 
to be the minimum projected separation from the center- 
of-mass, and normalize Ie to the total mass of the binary. 
For 6 > 1.0 and q < 10~^, we choose Mmin to be the min- 
imum separation from the more massive lens component, 
and normalize Ie to this component as well. For q — 1 and 
b — 1.5 we choose Umin to be the minimum projected sep- 
aration from the center-of-mass, and normalize Ie to the 
total mass of the binary. Finally, for g = 1 and b > 2.0, 
we choose Umin to be the minimum separation from one of 
the components, and normalize tE to the same component. 
These choices were made in order to approximate the tran- 
sition between close and wide binaries and minimize the 
deviation globally. Since this qualitative example does not 



use best-fit single lens parameters, S will be overestimated 
in all cases; we discuss this effect fully in § 5.1. To cover 
the full range of parameter space for which the detection 
efficiency is high, we choose 5 values of q, logarithmically 
spaced between 1.0 and 10~^, and b = 0.2, 0.6, 1.0, 1.5, 
2.0, and 3.0. The results are shown in Figure |[ where con- 
tours of i5 = ±0.01 and ±0.05 are plotted with the caustics 
{S = oo). Light curves are one-dimensional cuts through 
these diagrams. Since current survey programs rarely alert 
on-going events with Wmin- 1, only trajectories that pass 
within 9e (circles in Figure |^) of the primary lens will be 
observed by the monitoring teams. 

In Figure |[ light curves resulting from the sample tra- 
jectory in Figure]2|are displayed. This trajectory is rather 
typical, with Umin = 0.5 and a randomly chosen angle 6. 
While for q — 1.0 and 10""'^, most trajectories will exhibit 
considerable {S > 0.05) deviations, for qS 10~^, the mag- 
nitude of the deviation will depend strongly on the value 
of the 9. For example, although the light curve in Figure || 
for q — 10^"^ and b ~ 0.6 exhibits no significant deviation, 
inspection of Figure ^ reveals that a trajectory with the 
same impact parameter but 9 ^ 90° would be likely to 
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Fig. 3. — The fractional deviation, 5, from the single-lens light curve as function of time from the peak of the single lens light curve, to, in 
units of the Einstein ring crossing time t^, for the trajectories shown in Figure 2. 



have larger 6. This illustrates why a calculation of the 
detection efficiency of a non-anomalous light curve must 
involve an integration over 9, which is degenerate in the 
single-lens case. Figures ^ and || also illustrate why the 
detection efficiency of a light curve depends strongly on 
its itmin- The sample trajectory for q = 10~^ and b = 1.0, 
for example, does not deviate more than 6 ^ 0.05, yet al- 
most all trajectories with smaller Umin would exhibit much 
larger deviations. 

Several conclusions can be drawn from inspection of Fig- 
ures ^ and ||. First, for g ~ 1 and 0.2 < 5-3.0, nearly all 
trajectories have deviations S > 5%. Second, for all sep- 
arations 6, events with Umin~ 0.1 will have a much higher 
detection efficiency than larger Umin events. Third, for 
small mass ratios {qS 10^^), it is likely that only a small 
fraction of detected events will exhibit caustic crossings, 
since, for these mass ratios, the area covered by the caus- 
tics is considerably smaller than the area covered by the 
S = ±0.05 contours. Finally, for small q, detection effi- 
ciencies for light curves with typical impact parameters 
will be substantial only for companions with separations 
0.6- b~ 1.6 (i.e., in the 'Lensing Zone'). 

5. DESCRIPTION OF THE ALGORITHM 



An schematic of our basic algorithm is displayed in 
Figure ^. The top panel shows a simulated light curve, 
which for the moment we use as a stand-in for an ac- 
tual observed light curve. (For details on how this curve 
was generated, see § 5.1.) For simplicity, we assume 
that the baseline flux is known perfectly {Fq = 1) and 
that the event is not blended {fs — 0). We fit a re- 
stricted PSPL model with the three remaining free param- 
eters t_E,Uinin,io to this 'obscrvcd' light curve, and obtain 
xt = 0.96 for 300 degrees of freedom, with best-fit pa- 
rameters {tE — 15days, Mmin — 0.3, io — 15.0). This value 
of xt indicates that the observed light curve is consistent 
with the PSPL model. In order to calculate the efficiency 
function e(b,q) for this event, we must determine what 
fraction of all possible light curves arising from a b — q 
binary lens is incompatible with the observed light curve. 
The middle panel of Figure shows contours of constant 
S for a binary with q — 10 and b — 1.5, and a sample 
trajectory with Umin = 0.3 and angle 6 = 120° between 
the trajectory and the binary axis. Since this source path 
does not cross any regions of significant deviation, the cor- 
responding observed light curve would be consistent with 
PSPL within the precision of typical monitoring photom- 
etry. Light curves resulting from trajectories with other 
values of 9 would be inconsistent with the observed light 
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Fig. 4. — Top panel: The difference in magnitude, —2.5 logjQ(j4), as a function of time for a simulated microlensing event with photometric 
precision of 2% and with uniform sampling at the rate /200 from alert until ■ The dotted line shows the best-fit point-source point-lens 
(PSPL) event, which has parameters 4^; = 15 days,niiiin = 0.5, to = 15 days. Middle panel: Contours of constant fractional deviation S from 
the single-lens magnification, as a function of source position {rj, g), in units of the Einstein ring radius 6e- The solid line shows a trajectory 
with ttniin = 0.5 and 9 = 120°. The contours are the same as in Figure 2. Bottom panel: The difference in between the best-fit binary-lens 
Hght curve with b = 1.5 and q = lO'^ and the best-fit PSPL Ught curve, as a function of 6. The horizontal fine indicates a detection threshold 
of AXjjjj,j,|,j^ = 225. The detection efficiency e of the light curve in the top panel to companions with b = 1.5 and q = 10"'^ is the fraction of 
all possible trajectories for which Ax^ > ^Xthrcsh' '-'^ ^^^^ case, e = 5%. 



curve. In order to determine the incompatibility of a & — g 
binary lens model with the observed light curve as func- 
tion 0, for each fixed 9 we find the best-fit binary-lens light 
curve, but leaving UminjiB and as free parameters. We 
then calculate Ax^{b,q,0), the difference in between 
the best-fit binary-lens light curve and the best-fit single 
lens light curve as a function of 0. The bottom panel of 
Figure H shows Ax^ as a function of for a binary with 
q = 10 and b — 1.5. Since only a small fraction of all 
possible trajectories would give rise to binary-lens light 
curves that are statistically incompatible with the data, 
the detection efficiency of this light curve is small for this 
parameter combination. Quantitatively, the detection ef- 
ficiency is simply the fraction of all possible trajectories 
{0 < < 2tt), for which Ax'^{b,q,0) > AxLosh- Thus 
e{b, q) is given by, 



e{b, q)^— d0 Q[Ax\b, q, 0) - AxLcsh], (9) 



where Q[x\ is a step function. For the event depicted in 
Figure e = 5% for b = 1.5 and q = 0.001 and a detec- 
tion criterion of Axt^resh = ^^5. This process must then 
be repeated for all q and &, and then for all light curves 
that are consistent with the point-lens model. 

For Gaussian errors, x^ is the best measure of goodness- 
of-fit, and the significance of the detection can be altered 
by adjusting Ax^^^^^^, the minimum Ax^ between the 
best-fit single and binary lens models required for a de- 
tection. The choice of Ay^^^^^^ required for a 'detection' 
is arbitrary, but it should be kept in mind that error dis- 
tributions for actual monitored events are far from Gaus- 
sian, and usually contain systematic errors with unrec- 
ognized correlations at the few percent level. In light of 
this, we choose a rather conservative detection criterion, 
^Xthresh > 100. The choice of the appropriate detection 
criterion for realistic error distributions will likely depend 
sensitively on, and be determined by, the actual error dis- 
tributions themselves. We return to this point with a dis- 
cussion of the dependence of detection efficiencies on the 
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detection threshold in § 5.1. 

The basic steps in calculating the detection efficiency of 
events consistent with the PSPL to a specific b — q binary 
are summarized below. With the appropriate modifica- 
tions, a similar algorithm can be applied to the analysis of 
detection efficiency for any microlensing anomaly, includ- 
ing those due to binary sources, lens rotation and parallax 
effects. 

(1) Fit each event with a single lens model by minimizing 

(or some other suitable goodness-of-fit estimator). 
Evaluate for this model. 

(2) Hold the angular separation and mass ratio (5, q) 

fixed. For each source trajectory 9, find the binary 
lens model that best fits the observed light curve, 
leaving tE,tQ, and Wmin as free parameters. Evalu- 
ate the difference Ax^{b, q, 9) between the single- lens 
and binary-lens fits. 

(3) Find the fraction of all binary-lens fits for the given 

(5, q) that satisfy the detection criterion (e.g. Ax^ > 
^Xthresh)- "^^^^ IS the detection efficiency e{b,q) for 
this event for the assumed separation and mass ratio. 

(4) Repeat items (2) and (3) for all {b,q). This gives the 

detection efficiency for the ith event as a function of 
b and q, ei{b, q). 

(5) Repeat items (l)-(4) for all events. 

The steps itemized above assume that the baseline fiux 
Fq is known perfectly, and that the event is not blended, 
(/b =0). In reality, one must always fit for the base- 
line flux Fq and blend parameter Jb ■ If an event is truly 
blended, including Fq and Jb in the fitting procedures can 
have a strong effect on the computed detection efficiency of 
the resulting light curve. Similarly, the algorithmic outline 
above assumes that the source can be treated as point-like. 
Including finite source sizes can also have a significant ef- 
fect on the inferred detection efficiency of a given event. 
In order to obtain an accurate estimate of the detection 
efficiency, these effects must be included in the fitting pro- 
cedure, and either the light curve itself or other data used 
to constrain the blending and finite source size parame- 
ters. In order to clearly delineate the effects of blending, 
finite source size, and the choice of the detection crite- 
rion, we will first assume that the baseline flux is known 
perfectly, the blending is negligible, and the source can 
be approximated as a point-like. The effects of detection 
criterion, finite source size and blending on the detection 
efficiency e{b,q) are then explored separately in §§ 5.1, 6 
and 7, respectively. 

The detection efficiencies calculated in the prescribed 
way for non- anomalous events can be used in several ways: 

(1) to place quantitative constraints on the absence of 
planets of certain 6, q in non-anomalous lensing event; 

(2) to estimate the average detection efficiency e(6, q) for 
a given dataset; (3) to estimate e(6, q) for hypothetical 
datasets as a guide to future observational programs, and 
(4) as a proxy for the detection efficiency of observed 
anomalous events, for which additional challenges exist 
(see discussion in § 8). 



5.1. Application to Artificial Data 

In order to explore more fully the effect of the param- 
eters Minin, 6, q, and detection criterion on the detection 
efficiency, and to test the robustness of the algorithm, we 
generate artificial light curves and calculate their detection 
efficiency e(b, q). 

Each simulated event is assumed to be alerted a,t A = 
1.54 (the smallest amplificati on that the MACHO team 
will alert, Alcock ct al. 1997b| ), and then continuously ob- 
served at uniform intervals of ts/'^QO until either after 
the peak or 3tE after the peak. A more realistic light 
curve would contain gaps due to bad weather or other ob- 
serving conditions. In order to isolate as much as possible 
the intrinsic dependencies of the detection efficiency, we 
use uniform sampling, and also assume here /s = and 
Fq = 1. At each observation, a residual is drawn from a 
Gaussian distribution with a = cfqAq, where ctq = 1% or 
2%, and is given by Equation (0). These parameters 
are roughly consistent with the best sampling a nd pho- 
tometric ac curacy of the PLANET collaboration (Albrow 
et al. 1998). Hereafter, all results will be for light curves 



with (Jo = 2% and observations until Ie unless otherwise 
noted. Three values of the impact parameter, Umin — 0.5, 
0.3, and 0.1, are investigated. 

For each light curve, we calculate e{b, q) using steps (2- 
4) in § 5. When fi tting the binary, w e employ a downhill- 
simplex method ( |Press et al. 199'^ , which usually con- 
verges quickly and robustly to the minimum. The fitted 
single-lens parameters are used as an initial guess, since 
best-fit binary parameters are typically close to the single- 
lens values for small mass ratios q. 

Since the fitting procedure can be computationally ex- 
pensive, we sample 9 only at intervals of 27r/100; our effi- 
ciencies are thus limited to a resolution of Ae — 0.01. Al- 
though one would like to sample the (5, q) plane as densely 
as possible, we are again limited by computational ex- 
pense. Since planetary events are the primary interest 
of most monitoring collaborations, and nearly equal mass 
binaries (10~^ < 9 < 1-0) will have e ^ 100% in the Lens- 
ing Zone anyway, we choose to restrict our attention to 
q < 10^^. Furthermore, the effect of finite source sizes 
and blending, which we wish to investigate, will be sub- 
stantially less dramatic for q^ 10~^ than for lower mass- 
ratio systems. We choose three mass ratios: q — 10^^, 
10"'^ and 10~^. Using a different detection algorithm, Ben- 
nett & Rhie (1996) found detection probabilities of ~2% 
for q = 10~^. Since this is comparable to our resolution 
Ae — 0.01, we will not extend our analysis to mass ratios 
smaller than q = 10^^. We calculate e from b = 0.2 to 2.0 
at intervals of 0.2, and then again at 6 = 3.0. 

Previous explorations of planetary microlensing detec- 
tion probabilities have not used the best-fit binary lens 
light curve, which is computationally expensive to com- 
pute. Holding tE,tQ and Mmin fixed will cause one to over- 
estimate the detection efficiency in two ways. First, for 
large mass ratios q~l{)~^, the secondary cannot simply 
be treated as a perturbation to the primary light curve. 
The presence of the secondary will have a significant ef- 
fect on the global (averaged over all trajectories) values 
of the best-fit parameters. For example, a close (5 ^ 1) 
equal mass binary will have a time scale that is a factor 
of 2^/^ larger than an otherwise identical wide {b ^ 1) bi- 
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Fig. 5. — The detection efficiency, t{b,q) as a function of the instantaneous separation of the binary, b for three different mass ratios 
{q = 10~^,10~^, and 10~*) and impact parameters (itmin = 0.1,0.3, and 0.5). Three detection criteria are shown: Ax^ > 225 (thick solid 
line), Ax'^ > 100 (thin sohd line), and caustic crossing (dashed line). The efficiency e using Ax^ > 225, but without fitting the binary-lens 
light curve, is also shown (dotted line). Throughout it is assumed that the events are sampled uniformly with 2% photometry at the rate 
tB/200 from alert to ts- See § 5.1. 



nary. Second, for small mass ratios q~ 10^"^, the majority 
of the detection efficiency will arise from relatively small, 
(5-10%, deviations. These deviations can be suppressed 
below the detection crite rion if the fit is allowed t o adjust 
to compensate for them (Griest & Safizadeh 1998). In or- 
der to facilitate comparison between our results and previ- 
ous calculations, and to gauge the error induced by hold- 
ing the parameters ts, to, and Wmin fixed when calculating 
the difference between binary-lens and single-lens magnifi- 
cations, we have also calculated e(6, q) without finding the 
best-fit binary, assuming instead that tE, to and Umin are 
the same for the binary and single lens fits. For 6 < 1, we 
choose the center-of-mass as the origin of the binary; for 
5 > 1, we choose the position of the primary. 

In Figure ^ the detection efficiency e(5, q) is displayed 
as a function of the dimcnsionless angular separation b for 
events that are followed until t — Ie- Three different im- 
pact parameters (umin = 0.1,0.3, and 0.5) and mass ratios 
[q = 10^^, 10^'^, and 10"'*) are investigated using different 
detection criteria: Axthresh = 100' '^Xt^resh = 225, and 
the criterion that the trajectory must cross a caustic in or- 
der to be detected. Also shown is e(6, q) for Axthrcsh — ^25 



without fitting the binary light curve. For all events, we 
have assumed a photometric accuracy of 2% and that ob- 
servations are carried out from alert until He after the 
peak. 

Figure ^ illustrates several points. First, the difference 
between the detection efficiency calculated using the two 
different thresholds, Ax^^^y-^sh ~ 225 and 100, is small and 
approximately constant at Ae ~ 0.05 over most of the pa- 
rameter space considered. However, since the magnitude 
of e(6, q) decreases with decreasing mass ratio q, the frac- 
tional difference increases. The exact choice of Axthrcsh 
thus has little effect on e(6, q) for deviations well above 
the detection threshold. For perturbations near the detec- 
tion limit, however, such as those arising from companions 
with mass ratio q = 10~^, the detection efficiency can vary 
by a factor of two depending on the choice of detection cri- 
terion. Second, the error induced by not using the best-fit 
binary lens light curve can be substantial for q = 10~^, 
because such companions cause significant anomalies over 
a large fraction of the light curve so that the parameters 
of the best fit binary- and single-lens models can differ 
dramatically. This is especially true for high magnifica- 
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Table 1 Lensing Zone Detection Efficiencies eLz('z): Point Source 



tion events with b > 1, for which the fractional deviation 
S from the best single-lens model depends critically on the 
choice of the binary origin, which can be significantly dif- 
ferent from the primary lens position depending on the 
value of b. The error induced is also large for q = 10"**, 
because these deviations are very near the detection limit. 
For q — 10~^, the deviations caused by the companion can 
be treated as a perturbation to the primary light curve, 
and are well above the detection limit, minimizing the ef- 
fect of not using the best-fit binary; we caution, however, 
that this is unlikely to be true for all realizations of re- 
alistic sampling. We conclude that one can avoid fitting 
the binary lens light curve for mass ratios q~ 10~^ only 
if the deviations are well above the detection limit of the 
observed light curve. Finally, as noted by Gould & Loeb 
(1992), we find that caustic crossing events are likely to 
comprise only a small fraction of all detected events. This 
is especially important because non-caustic crossing events 
are more prone to degeneracies and thus the most d ifficult 
to characterize ( Glaudi fc Gould 1997 ; |Gaudi 1998 ) . 

In order to quantify the effects of the various detection 
criteria and impact parameters on e{b,q), we tabulate in 
Table average detection efficiencies eLz('7) for the curves 
in Figure ^ integrated over the Lensing Zone (where the 
detection efficiency is the highest), 0.6 < b < 1.6, 



eLz(g) 



1.6 



dbeib,q) 



(10) 



The accuracy of these results are limited by the fact that 
we sample e(6, q) only at intervals of Ab = 0.2 in this zone. 
For any given event, however, the results are more se- 
cure, so that comparisons between elz for (umin = 0.1, g ~ 
10~^) and (umin — 0.1, g — 10~^) should be more reliable 
than comparisons between (umin — 0.1,(7 = 10^^) and 
(umin = 0.3, g = 10~^). Table |] illustrates that the frac- 
tional error induced by not fitting the binary light curve is 
smallest for q = 10"^, < 15%. For q = 10"^ and q = IQ-^ 
however, the error can be considerably larger, ^ 20% for 
Umin — 0.3 and q — 10^^, and ^ 30% for w,„i,i — 0.5 
and q = 10^*. The fractional difference in elz between 
'^Xthrosh — 100 and 225 can be substantial, especially for 
q = 10"'*, where it is always greater than 50%. Finally, 
caustic crossing anomalies comprise a relatively small frac- 
tion of the total events, representing at most 35% of the 
integrated detection efficiency, and decrease in importance 
for large impact parameters and smaller mass ratio. 



If all lensing primaries have planets distributed uni- 
formly in b, the numbers in Table ^ represent the fraction 
of all events with the given Umin that would exhibit de- 
tectable deviations with the given q and 0.6 < b < 1.6. 
For larger mass ratio companions with q — 10~^, the de- 
tectable fraction is quite large, eLz(g) ~ 30%, and remains 
substantial, eLz(g)~10%, even for "Jovian" companions 
with q = 10~^. For small companions with q = 10"'^, 
however, the detectable lensing zone fraction drops signifi- 
cantly below 10% in all cases, although the exact numbers 
are somewhat uncertain due to the poor sampling in b. 
We conclude that microlensing will be able to place strong 
constraints on the frequency of double lenses with q~ 10"^ 
and mild constraints on systems with q = 10~^ compan- 
ions, but will be unable to meaningfully constrain systems 
with q~ 10""*, unless the sampling and photometric preci- 
sion are significantly better than those assumed in these 
simulations. 

6. FINITE SOURCE EFFECTS 

The calculations in § 5 implicitly assumed that the mi- 
crolensed source was point-like, so that the magnification 
of the source is infinite at the caustics. The magnification 
Afs of a source with finite size is given by the integral of 
the point-source magnification over the face of the star. 



AUt) 



Jd'^rA{t;r)I{r) 
/d2r/(r) 



(11) 



and is equivalent to the intensity-weighted area of the im- 
ages (numerator) divided by the intensity-weighted area 
of the unlensed source (denominator). Here /(r) is the in- 
tensity profile of the source. The finite size of the source 
smoothes and broadens the discontinuous jumps in mag- 
nification near caustics by an amount that depends on the 
angular size of the source, 0,, in units of 9e, 
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(12) 

where i?* is the physical size of the star. For the scal- 
ing relation on the right of equation (|l2[), a lens distance 
of DoL = 6 kpc and a source distance of Dqs = 8 kpc 
has been assumed. For a given , uniform sources 
(/ = constant) will have a larger effect on the magnifi- 
cation than limb-darkened sources. Since we are inter- 
ested primarily in the magnitude of the effect that finite 



13 



1 1 1 1 1 1 1 1 

h — n H 


1 1 1 1 1 1 1 1 1 1 1 










/ b-1.2 
- ~/ p* = 

1 /■ 1 1 1 1 1 1 1 1 1 


® p. = 0,03 - 



-0.5 0.5 



Fig. 6. — Contours of constant fractional deviation, S, from the single-lens magnification for a point-source (left panels) and source with 
radius p* = 0.03 (right panels), as a function of source position (»;, ^), in units of the Einstein ring radius, 9^, for mass ratio q = 10~* and 
dimensionless separations b = 0.8, 1.0, 1.2. The contours are the same as in Figure 2. The light curves arising from the sample trajectories 
(solid and dotted straight lines) are shown in Figure 7. 



sources will have on the detection efSciency, we will as- 
sume a (less-realistic) uniform source profile, which also 
increases computational speed. 

Since the caustic of a point lens is a single point at the 
center of the lens {u = 0), the magnification of a finite 
source will differ from that of a point source only when 
the source approaches the center of the lens, Umin ^ P*- 
For typical sources and lenses, -C 1 (c.f. equation (p^), 
so that for single lenses finite source effects are noticeable 
only in high magnification events (umin 1), for which 



A 



(equation 



Assuming a uniform source, the 



single lens magnification can be found analytically in this 
limit ( ISchncidcr, Ehlcrs, fc Falco 19921 ), 

(13) 



where 



Biz) 



AisAt) = Aoit)B[z{t)], 



^E{z) z < 1 

lz[Eil/z)~il^z-^)Kil/z)] z>l 



(14) 



Here K and E are the complete elliptic integrals of the 
first and second kind, and z = u/p^. 

For a binary lens, the finite size of the source affects the 
magnification whenever the source approaches a caustic 



or, more precisely, wherever the second derivative of the 
magnification is large. As can be seen in Figure |2[ for bina- 
ries with q ^ 1 much of the region inside the Einstein ring 
satisfies this condition: caustic approaches and crossings 
will be common. On the other hand, the results of § 5.1 in- 
dicate that caustic crossings make up only a small fraction 
of all detectable events for binaries with mass ratios con- 
sistent with planetary systems {q~ 10~^). Nevertheless, as 
can be seen from Figure ||, in order to produce a detectable 
event trajectories must pass close to caustics, where the 
gradient of the magnification is large and finite source ef- 
fects are non-negligible. Finite source magnifications for 
binary lenses cannot be found analytically. Numerical in- 
tegration of the point source magnification over the face 
of the star is difficult, as the divergent magnification near 
caustics causes the results to depend critically on the in- 
tegration grid size. A more robust method is to compute 
the total area of all images and then divide by this area 
by that of the source to find the magnification of the fi- 
nite source at that position. Numerous methods have been 
sugge sted: we will integrate over the boundary of the im- 
ages (Kayser &: Schramm 1988; Gould fc Gaucherel 1997; 
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Fig. 7. — The fractional deviation, S, from the single-lens light curve as function of time from the peak of the single lens light curve, to, 
in units of the Einstein ring crossing time tg, for the trajectories shown in Figure 7. The solid line is for a point-source, and the dotted lines 
are for a source of radius p« = 0.03. 



Dominik 1998). For alternative methods, see Bennett & 
Rhie (1996), Wambsganss (1997), and Griest & Safizadeh 
(1998). 

What is relevant to this discussion is not the difference 
between the finite source and point source magnification, 
but the effects of finite source size on the determination of 
the detection efficiency e. For nearly equal-mass binaries 
{q~ 10^"'^), the magnification may be altered considerably 
by finite source effects without substantially altering e. 
This can be seen by comparing the size of the deviations 
(6) in the q = 1.0 and q = 10~^ panels of Figure ^ to the 
size of a large source (p* = 0.03). For mass ratios as small 
as q~ 10~^, however, the size of the 6 structures is compa- 
rable to that of a large source, and finite source effects are 
important to a proper determination of the detection ef- 
ficiency e. Roughly speaking, finite source effects become 
important whenever the source size becomes comparable 
to the Einstein ring of the companion, 9^,~0p — q^l'^dEi 
or, 

1/2 



(15) 



This criterion is satisfied at ~ 0.1 for q ~ 10~^, whereas 
for q ^ 10"'^ source sizes ~ 0.03 will begin to seriously 
affect e. Since the largest sources routinely monitored in 



the Galactic bulge are clump giants, with p ~ 0.03, finite 
source effects will be negligible for q > 10^"^, but must be 
considered for smaller mass ratios. 

Although the magnitude of the perturbation will always 
be suppressed in the presence of a finite source, it will 
also be broadened. Finite sources thus have competing ef- 
fects on the detection efSciency e: e is decreased because 
previously significant deviations are suppressed below the 
detection threshold, but increased for those trajectories 
for which the limb of the star grazes a caustic (or high 
magnification area) yielding a significant deviation where 
no significant deviation would have occurred for a point 
source. The net result of these two competing effects will 
depend on the specific value of q and h. 

Contours of constant fractional deviation 5 of a binary 
from a single lens magnification, are illustrated in Fig- 
ure ^ for both a point source and a finite source. We 
choose — 0.03, a relatively large source, and q — 10"'', 
the smallest mass ratio we consider, in order to present 
a scenario in which the source size will have an extreme 
effect. For the finite-source cases, the fractional devia- 
tion is computed with respect to a finite-source point-lens 
magnification, as given by equation ( p!3| ) . Figure || clearly 
demonstrates that, for this mass ratio and source size, the 
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Fig. 8. — The detection efficiency e as a function of the dimensionless separation of the binary, b, for three different mass ratios q and 
three events with different minimum impact parameter Mmin- The sohd hne is e for a point-source and a detection criterion Ax^ > 225; the 
dotted line is e using the same detection criterion for a source of radius p* = 0.03. It is assumed that the events are sampled uniformly with 
2% photometry at the rate t£;/200 from alert to t^- 



differences between the point and finite source magnifi- 
cation are dramatic. Both the shape and size of the S- 
contours are altered considerably. 

As can be seen from Figure 0, planetary perturbations 
with 6 < 1 are qualitatively different than those with 6 > 1 
(Gould & Loeb 1992; Gaudi & Gould 1997; Wambsganss 
1997). Consider the case b — 0.8: the perturbation is sub- 
stantially depressed by the finite source and the S = 5% 
contours have nearly disappeared. This is because, for 

5 < 1, regions of constant positive and negative deviation 
are closely spaced and of nearly equal area so that the 
smoothing induced by a large source tends to cau se a can- 
cellat ion leaving a deviation that is n early zero (Bennett 

6 Rh ^e 1996| ; pould fc Gauchcrcl 1997| ). The effect is even 
more prominent for 6^1, where the regions of positive 
and negative deviation are especially closely spaced (Ben- 
nett & Rhie 1996). Obviously, for these two parameter 
combinations, the planet is unlikely to be detected and 
e — > 0. For b > 1, regions of positive deviation encompass 
considerably more area than those of negative deviation 
(at a fixed value of |(5|~5%), and the cancellation is less 
dramatic. As a result, e(6 > 1) will be less affected by 
finite sources than e(6 < 1). In both cases, the perturba- 



tions caused by the central caustic (near ( = 0,7] = 0) have 
dropped below S — 1%. Central caustics are an important 
channel to planet detection in high magnification events; 
detection e fficiencies for these event s will be highly sensi- 
tive to p.^, ( priest fc Safizadeh 199^ ). Fortunately, this is 
a class of events for which can often be measured. 

In Figure]^, the light curves resulting from the trajecto- 
ries shown in Figure ra are displayed. The trajectories were 
chosen to create a significant point source fractional devia- 
tion, but are otherwise representative. Dramatic cancella- 
tion can be seen in the light curves for b = 0.8 and b = 1.0. 
For b = 1.2, the deviation is substantially suppressed but 
is also broader. For photometry of sufficient precision, the 
detection efficiency e for this parameter combination will 
actually be increased. 
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Table 2 Lensing Zone Detection Efficiencies eLz('z), 
Cthresh " 225: Finite Source 



Ax? 



To malce a quantitative comparison, we have calculated 
e in the same manner and for the same parameters as in 
§5.1, but now compare the simulated light curves to the 
best-fit finite source binary light curves with = 0.03. 
The results with a detection criterion Ax?jjj,^gj^ = 225 are 
shown in Figure ^ along with the corresponding point 
source efficiencies from Figure |[ In agreement with the 
estimate from equation (15|), the detection efficiency e for 
q = 10~^ companions is nardly affected. For q = 10"'^ 
and separations 6 < 1, e is similar to or smaller than finite 
source efficiencies, but for wider separations & > 1, e can 
be either somewhat smaller or somewhat larger due to the 
finite source size. The difference is dramatic for q = 10^'*: 
source sizes corresponding to bulge giants always yield ef- 
ficiencies 5 1%. As in § 5.1, we calculate lensing zone effi- 
ciencies elz (equation (p^); the results are shown in Ta- 
ble ^ We conclude that finite source sizes have negligible 
effect on e for 10~^, but sources as large as bulge giants 
(p* ~ 0.03) can have a dramatic effect for smaller com- 
panions, either increasing or decreasing e (g ~ 10"'^), or 
wiping out the detection efficiency completely {q~ 10^^). 

Unfortunately, for individual events, the value of is 
very poorly constrained. While it is possible to estimate 
the physical size of the source from its color and magni- 
tude, this cannot be translated to the dimensionless pro- 
jected size if the value of 9e remains unknown. The 
detection efficiency e for most events could be in error 
therefore by many tens of percent (c.f. Figure |^ and Ta- 
ble ^). We discuss method of dealing with this difficulty 
in 5 8^. 



7. BLENDING 

A microlensing event is blended whenever unresolved, 
unlensed background light contributes significantly to the 
observed baseline fiux of the source star (i.e. Jb ^ 0). 
From equation (|l|) , the observed magnification in the pres- 
ence of blending, As, is related to the true magnification 

by, 

AB^{l + fB)-\A + fB), (16) 

where Jb is the ratio of the blend flux to the true lensed 
source fiux. 

Blending will have two effects on the detection efficiency. 
The first is a suppression of the deviation caused by the 
binary. From equation (|6|) and equation (|l|) , it is straight- 
forward to show that the fractional deviation in the pres- 
ence of a blend, i5b, is related to the true fractional devi- 



As before, Aq is the magnification of the best-fit PSPL 
light curve. Note that equation (|l^) applies to any 
anomaly that produces a deviation from the standard 
PSPL light curve, including parallax, binary source, and 
finite source effects. Since is a function of time, the 
magnitude of the suppression will also be a function of 
time, such that deviations occurring closer to the peak 
{u ~ Mmin) wiU bc less suppressed than those occurring 
near the beginning and end of the event. Figure |^ shows 
the ratio 6b /S as a function of the angular separation u of 
the lens and source, for three values of Jb corresponding 
to relatively mild blending, /s < 1. Anomalies occur- 
ring near the peak of high ma gnification events (u ^ ) 
will be only slightly suppressed (Griest & Safizadeh 1998), 
while repeating events caused by wide binaries (Di Ste- 
fano & Mao 1996; Di Stefano & Scalzo 1999b) can be sup- 
pressed by as much as (1 -I- /b)~^. Overall, the amplitude 
of the suppression for mild blending is relatively small, 
Sb/S < (1-1- /b)~-^, and thus will not have a large effect 
on the detection efficiencies for binaries with mass ratios 
q > 10^^ for which the fractional deviations are usually 
large. For mass ratios consistent with planets, q < 10~^, 
however, a substantial fraction of detected events will have 
maximum fractional deviation (5i„ax~0.1 so that even a 
suppression of 0.5 can have a significant effect on plan- 
etary detection efficiencies. 

The second effect that blending has on detection effi- 
ciencies is to alter the presumed distribution of Umin in 
observed events. The intrinsic distribution of Umin is flat 
up to the magniflcation threshold set for detection by the 
survey teams (e.g., A > 1.34 corresponds to Umin < !)• 
This threshold is calculated in real time near the begin- 
ning of the event when a robust determination of blending 
is not possible. Consequently, the magniflcation of the 
event at any time is assumed to be the total flux divided 
by the basehne flux, Ab — {I + fB)~^{A + Jb), which is 
less than the true magniflcation A for all non-zero blend- 
ing values. Thus, in the presence of blending, an event 
will require a larger magnification and thus smaller intrin- 
sic Umin in order to pass the detection criterion, which is 
per force ap plied to the observed quantity Ab (see Alcock 



et al. 1997a for an example and discussion). Since the in- 
trinsic detection efficiency is larger for smaller impact pa- 
rameters (see § 5.1 and Figure ||), the blended event will 
be more sensitive to the presence of planets than would be 
calculated for a Umin based on the (erroneous) assumption 
that /s = 0. Blending thus affects binary detection effi- 
ciency in two competing ways: the efficiency is decreased 
by suppressing the amplitude of observed deviations, while 
at the same time it is increased due to the skewing of the 
observed Umin distribution to smaller values. The net ef- 
fect on the detection efficiencies will depend on the values 
of b and q and vary on an event-by-event basis. 

Clearly, blending must be considered when calculating 
binary and planetary detection efficiencies. Since blending 
effects are relatively easy to quantify, this poses no seri- 
ous complication as long as the blending parameter, /b, 
can be accurately determined for individual events. Un- 
fortunately, as discussed by Wozniak & Paczynksi (1997), 
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Fig. 9. — Ratio of blended Sg to true fractional deviation 5, as a function of impact parameter u for three different blend fractions, fs- 
Microlensing events are generally alerted when m < 1. 



blending can be extremely difficult to determine for indi- 
vidual observed light curves, due to the serious correla- 
tions in the parameters Fo,tE,Umin and fs in the pres- 
ence of blending. The degeneracy is especially severe in 
two regimes. In the 'spike' regime, defined by /s ^ 1 and 
Wmin <SC 1, ifi and Umin caunot be measured separately and 
only the deg enerate combination i£;/itmin is measurable 
(Gould 1996). This regime is most important in severely 
crowded fields, such as those towards M31. The second 
regime, defined by large impact parameter (wmin ~ 0.5) 
events with modest blending /s < 1, is more common in 
bulge fields, and thus the focus of our current attention. 

In order to illustrate the difficulty in quantifying the 
fraction of blended light, we calculate the range of allowed 
values inferred for fs from fits to our fiducial simulated 
light curves (i.e., observations from alert until ts, 2% er- 
rors, and no blending). Fixing fs at some value, we find 
the best fit to a point-lens model allowing the parame- 
ters Mmin, i_E, ^Oj ^ud _Fo to vary. We then compute Ax^ 
between this fit and that assuming /b = 0. The result- 
ing Ax^ as a function of fs is shown in Figure |l^. For 
itmin = 0.1 and Uinin = 0.3, the degeneracy is not severe; 
the 3cr allowed ranges in Jb are Jb < 0.125 and fs < 0.45, 
respectively. For Umin = 0.5, the blend fraction is almost 
completely unconstrained because for large Mmin events 
constraints on blending arise mostly from the combination 
of information from the wings (|t — ^ol ~ ^b) of the event 
and the baseline (|t — io| ^ i^). Thus, without a baseline 
measurement, the value of Fq can be arbitrarily adjusted 
to compensate for large values of /s without significantly 
affecting the fit. With this in mind, we have also com- 



puted the same Ax^ statistic for a simulated light curve 
with observations from alert until 3tE- Here the blending 
is much better constrained, fs < 0.27, with additional im- 
provement if the errors are reduced by half, in which case 
Ib < 0.15. Since the majority of the constraint comes 
from sampling the wings and baseline of the light curve, 
it would be more efficient for monitoring teams to concen- 
trate on more, rather than better, measurements. 

The right panels of Figure ^ demonstrate how an 
inaccurately-determined blend fraction can affect the de- 
termination of Umin and thus the detection efficiency e. 
Here we show the ratio between the value of Wmin deter- 
mined by assuming a constant blend fraction /s, Umin.fj 
and the true value Umin- The value of Umin deduced for 
blend fractions between /s = to /s = 0.4 can vary 
by nearly 50% resulting in quite different inferred detec- 
tion efficiencies (c.f. Figure ph. To quantify this, we have 
calculated e as a function of /b for an extreme example, 
Umin = 0.5 and observations from alert until Ue- The 
procedure for calculating e is the same, except now /b is 
fixed at an assumed value and Fq is included as a free pa- 
rameter in both the single lens and binary lens fits. The 
results are shown in Figure |l^, where we plot eiq = 10~^) 
as a function of b for = (same as Figure |p, 0.2, 0.4, 
and 0.6. Recall that all these blend fractions are statisti- 
cally indistinguishable for this light curve. The differences 
in e arc dramatic. The suppression of binary anomalies 
induced by blending of /s = and 0.2 causes a drop in 
detection efficiency for separations 6-0.6 and b~ 1.6. In- 
side the Lensing Zone, however, the net effect is a dramatic 
increase in e due to the lower value of Umin required to pro- 
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Fig. 10. — Left panels: The difference in between a point-lens fit with no blending and one with the blend fraction held fixed at fg, for 
observations from alert {A = 1.54) until t = (dotted line) and from alert until t = 3tE (solid line). The effect of improving the photometric 
precision by a factor two while monitoring to t = Stg is also shown (long-dashed line). Right panels: The ratio of the fitted value of u„^i„ to 
the true value of Umin as a function of the blend fraction. Line types are the same as the left panels; the long-dashed line is coincident with 
the solid line. 



duce the observed light curve for increasing /g. This can 
be appreciated best by examination of Table 0, which tab- 
ulates clz as a function of fs- To the extent that these 
values of fs cannot be distinguished from one another by 
the light curve alone, they are all equally likely, and thus 
e can be quite uncertain. 
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Table 3 Lensing Zone Detection Efficiencies eLz(9 — 
10^2). Blending 

Blending can give rise to a serious uncertainty in de- 
rived detection efficiencies but, unlike finite source size, 



the blend fraction can be determined with sufficient accu- 
racy for most events. To do so, however, requires precise 
measurements by the monitoring teams during the wings 
of the event and at baseline. Without a reasonable quan- 
tification of the blend fraction, the detection efficiency of 
individual events will be very uncertain. 

8. APPLICATION TO REAL DATA 

In previous sections, we used artificial data to explore 
several effects that can influence significantly the deter- 
mination of the detection efficiency e of individual light 
curves, including detection criteria (§ 5.1), finite source 
size (§ 6), and blending (§ 7). These effects are often dif- 
ficult to quantify in real data, for which sampling and 
photometric precision is likely to vary on an event-by- 
event basis, and with observing conditions and microlens- 
ing phase for individual events. Furthermore, real-time 
observational decisions may alter (increase) the sampling 
of clearly anomalous events from that of apparently non- 
anomalous PSPL events. 

8.1. Variable Sampling and Photometric Precision 
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Fig. 11. — The detection efficiency e as a function of the dimensionless separation of the binary, b, for mass ratio q = 10"'^ and four 
different blend fractions fs for an event with 2% photometry, uniform sampUng from alert until at the rate t£;/200, and minimum impact 
parameter u^^^ = 0.5. 



Since the algorithm we presented in § 5 uses the actual 
light curve, and thus the actual sampling and photometric 
uncertainties associated with each observed event, irregu- 



lar sa npling and variable precision arc taken into account 



on individual data points are often taken to be the formal 
errors reported by the PSF-fitting alg orithms of photo- 
metric reduct ion packages like DoPhot OSchcchtcr, Mateo 



explicitly in the determination of the detection efficiency e. 
The efficiencies based on artificial data presented in pre- 
vious sections to illustrate general principles will not be 
strictly applicable to real microlensing events, for which 
weather and other considerations prevent continuous mon- 
itoring with a sampling of i£;/200. In general, the effect of 
reduced or incomplete sampling will be to lower detection 
efficiencies. For extremely long events, or those (including 
very short) events alerted post-peak by the survey teams, 
substantial portions of the light curve will have no (dense) 
monitoring at all, and the effects on the detection efficiency 
can be quite devastating. In such partially-monitored light 
curves, the PSPL fit parameters {umin,tE,to, Fq and /b) 
will be very uncertain, and in extreme cases almost com- 
pletely unconstrained. Since the detection efficiency de- 
pends on these parameters, the resulting uncertainty in e 
will be quite large. Unless additional information is avail- 
able (e.g. from the survey teams) to constrain the fit, 
these data will add almost nothing to our knowledge of 
the abundance of planets since their detection efficiency 
cannot be reliably quantified. 

For Gaussian, uncorrelated errors, the statistic can 
be used as a measure of goodness-of-fit. Real measure- 
ment uncertainties are seldom truly Gaussian, especially 
in crowded microlensing fields where systematic effects as- 
sociated with seeing, scattered light, and detector char- 
acteristics become increasingly important. Uncertainties 



fc Saha 1993), which of ten underestimate the true scat- 
ter (Albrow et al. 1 



, Image subtraction te chniques 

( |Tomaney fc Grotts 199^ ; |Alard fc Lupton 1998| ) may al- 
leviate some of these difficulties, but for the moment are 
too cumbersome and slow to implement for multi-site, real- 
time reduction of large fields and thus have not yet been 
implemented by monitoring teams. An empirical correc- 
tion to account for the correlation of measured photomet- 
ric magnitude with the FWHM of the point spread func- 
tion often results in a more Gaussian error distribution 
whose average magnitude corresponds more closely to the 
form al DoPhot-reported error (Naber, private communica- 
tion; Albrow et al. 199£ ). As long as the detection criterion 
^Xthrosh maintained at a suitably high value (§ 5), the 
exact error distribution is likely to have little effect on the 
computed efficiency e. Remaining doubts can be assuaged 
by attaching the observed error distribution derived from 
constant stars to artificially generated PSPL light curves 
to calibrate both the 'false alarm rate' and efficiency of 
true detections with a given 



thresh 



criterion. 



8.2. Finite Source Effects 

Finite source effects pose a significant challenge to the 
robust determination of the detection efficiency e because 
the dimensionless source size, is unknown a priori. As 
we have shown in § 5, the finite size of the source should 
have negligible effect on e for g~ 10"^ but can have a sig- 
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nificant effect for q-lO^"^. Thus, without additional in- 
formation about , e can be determined robustly only for 
10~^, which is unsatisfactory for microlensing monitor- 
ing programs whose primary goal is to learn about small 
mass-ratio systems. A first-order estimate for could be 
obtained by measuring the angular size of the source, 0*, 
and a ssuming that the relative proper motion of the lens, 
II = V l/£'ol7 is equal to the mean relative proper motion 



of these two i mages traces out an ellipse who size is 9e 
( Walker 1995| ), thus requiring an an astrometric accuracy 
considerably smaller than 1 mas. Preliminary studies have 
shown that SIM, with its planned 4/ias accuracy, should 
be able to measure 6e reliably for a lmost all known mi- 
crolenses toward the galactic bulge (Boden ct al. 1998a 
iJeoii^, Han & Park 1999|; pominik & Sahu 19991; IGould 



(/i) for all lenses toward the bulge. The dimensionless size 
of the source would then be given by 



Salim 1999). Ground-based interferometers currently be- 



P* 



(18) 



The angular size of the source can be estimated by its 
color and magnitude, or by obtaining a precise spectral 
type through more resource-intensive spectroscopy. Un- 
fortunately, the value of (/i) depends on the assumed ve- 
locity and spatial distribution of the lenses. Moreover, 
even within a given model, the distribution of /i is wide, 
having a variance of a fac tor of ~ 3 and long tails toward 
higher and lower values (Han & Gould 1995). The true 
value of n thus could differ substantially from (fi) , leading 
to a large error in e for individual events. When aver- 
aged over many events, these errors should approximately 
cancel, but current monitoring programs are very far from 
this regime, especially for anomalous events. 

A somewhat better estimate for the effect of on de- 
tection efficiencies could be made as follows. Assume that 
an event has measured time scale Ie and angular source 
size 0* . For an assumed model of lens distances and veloc- 
ities, the expected distribution of proper motions, Qi{p), 
can be computed. The individual detection efficiency can 
then be approximated as. 



(19) 



where — 6^/ptE, Pmax is the maximum proper motion 
allowed by the observed light curve, and Pmin is some rea- 
sonable lower limit. This model-dependent estimate of the 
detection efficiency should be more accurate whenever fi- 
nite source effects are large, but it is time consuming to 
compute because e(6, q) must be determined for many p*. 

Ideally, one would like to determine p* directly for each 
individual lens. This can be done by measuring p, from 
singl e-color light curves for only 5% of events ( Gould] 



19941 Witt 1995), namely those with high peak magnihca- 
tion. With both optical and infrared photometry, p^ could 
be determined for app r oximately twice as many events 



dGould fc Welch 1996| ; [Witt 1995| ). If the lens is lumi- 
nous, one could measure the proper motion of the lens 
directly using accurate astr ometry and a high res olution 
instrument, such as HST ( paudi fc Gould 1997| ). This 
will not be possible for all events, however, and requires a 
long temporal baseline. The angular Einstein ring radius 
9e can be determined directly by measuring the centroid 
shift of the two unresolved images created by the lens. As 
the lens passes across the line of sight to the source, these 
two images move and change magnification. The centroid 



ing developed, such as the Keck Testbed Interferometer, 
should be able to measure 9e for a smaU er, but substantial, 
fraction of events (Boden et al. 1998a). Measurements of 
this kind would require coordination and cooperation be- 
tween microlensing and astrometric communities, but the 
results would be well worth the effort. 

In summary, if microlensing monitoring teams can make 
use of other resources, such as HST and SIM, the best 
method of dealing with finite source size effects is simply 
to measure the dimensionless source size p, for each indi- 
vidual event directly. In the absence of these options, the 
effect of finite source size on e must rely on statistical and 
model-dependent estimates of the distribution of p,. 



8.3. Blending 

As discussed in § 7, the accurate determination of the 
blending, Jb, in individual light curves is essential to the 
accurate determination of the detection efficiency. Precise 
measurements during the wings of the event and at base- 
line can allow the quantification of fs , but may not always 
be possible: bad weather may prevent wing measurements 
in some events and the faintness of the source star may 
preclude precise baseline measurements in others. Im- 
proved or alternate methods of quantifying blending would 
be beneficial; several have already been suggested. 

Since the blending problem arises mainly because the 
PSFs of individual stars are blended together in crowded 
microlensing fields, simply improving the spatial resolu- 
tion of the observations eliminates or reduces the blended 
light in most cases. The resolution needed is roughly that 
of HST (Han 1997), for which only a modest use of re- 
sources is required. Nevertheless, even the resolution of 
HST will be insufficient to resolve any blended light that 
might arise from the lens itself or unresolved companions 
to the lens or source. Blending caused either by the lens 
or unrelated projected stars near the source or lens will 
generally cause a color shift in the combined light as the 



source is magnified during the course of the event (BuchaL 
ter & Kamionkowski 1996). Since most stars in the bulge 
have nearly the same color, the color shift is expected to 
be small and difficult to measure. Alternatively, the blend 
fraction can be quantified by measuring the centroid shift 



of the blended PSF during the co urse of the event ^ ( Gold- ^ 
berg & Wozniak; poldberg 1998 ). Blend fractions can be 
reliably determined in this way only for heavily blend ed, 
high-magnification events (Han, Jeong, & Kim 1998) in 
which the centroid shift is prominent. Finally, improve- 
ment can be made through using image subtraction rather 
than the usual PSF-fitting photometry. Image subtraction 
reveals only the time- variable flux, AF = Fo[A{t) — 1], so 



^This is not to be confused with the 0(1 mas) centroid shift caused by the motion and variable magnification of the two images created 
during a microlensing event. The centroid shift due to blending is ~ d{AB — 1)/Ab, where d is the separation of the blended sources and Ab 
is the observed (blended) magnification. Since d is of order the resolution, 1", this centroid shift is measurable using current ground-based 
observations. 
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that constant sources of blended light are removed. Al- 
though image subtraction does not remove the correlations 
between fitted parameters that make blending problem- 
atic, it is useful since photometri c uncertainties are gener- 



ally near the photon n oise limit (Tomaney & Crotts 1996 



Alard & Lupton 1995), which allows better discrimination 



between subtle differences in blended light curves. 
8.4. Anomalous Events 

If solar systems like our own are not atypical, dense, 
precise photometric monitoring will eventually result in 
the detection of anomalous events consistent with plane- 
tary perturbations. The algorithm of § 5 can be used to 
place direct limits on planetary systems in non-anomalous 
events and give meaning to the large number of null results 
now in hand, but what of the detection efficiency of events 
in which planetary anomalies are actually observed? 

By the detection efficiency ei{b,q) of an anomalous 
event, we mean, by analogy with the definition for non- 
anomalous events given in § 4, the probability that a lens- 
ing companion with parameters b, q would be detected in 
light curve i (using the same Ax^jjj,^gjj criterion) after inte- 
grating over source trajectory angles 9. In principle, if the 
parameters to,tE,Uniin, Fq and fs can be well-determined 
in a binary fit that meets the Axthrosh criterion, e(6, q) can 
be calculated in a manner similar to that presented in § 
5. A PSPL light curve with these parameters can be gen- 
erated with the actual sampling and photometric errors 
of the observed event and the method of § 5 applied for 
the b, q combination that produced the best-fit 'detected' 
binary. 

Several difficulties inherent to anomalous light curves 
must be addressed. First, for large mass ratio {q~ 10"^), 
the source trajectory 9 may actually be constrained by 
the light curve itself, but in a way that is partially de- 
generate with other fitting parameters. Second, even 
for small mass ratios, i.e., true planetary anomalies, 
the parameters to,tE,u^in, Fq and fs may not be well- 
constrained, though of course this difficulty is present for 
non-anomalous events as well. Third, the anomaly, if large, 
may produce changes in magnification significant enough 
to alter the photometric precision obtained at the phase 
of the perturbation. Most planetary anomalies will have 
rather gentle changes in magnification (Figure |) so that 
this effect may not be severe for g< 10^"^. Fourth, if de- 
tected real time, observer intervention may alter the pho- 
tometric sampling in a way that is also phase-dependent 
(i.e. higher sampling after the anomaly detection than be- 
fore). These last two possibilities present difficulties when 
integrating the corresponding PSPL light curve over source 
trajectory 9, since the altered sampling and photometric 
errors will not correspond to the phase of actual anomalies 
for most choices of 9. 

One approach to handling this increased sampling at the 
phase of the anomaly is to resample the 'post-anomaly' 
portion of the light curve in a manner that is consistent 
with what would have been the sampling had the anomaly 
gone unnoticed. (Determining how to do this in the pres- 
ence of weather and other observing facts-of-life may not 
be trivial) . Both the initial test of binarity and the proce- 
dure to determine the detection efficiency would then be 
performed with the sparser, resamplcd light curve. The 
full dataset would be used only to refine the best-fit pa- 



rameters. Alternatively, a Bayesian approach to analyzing 
the full dataset (anomalous and non-anomalous) could be 
employed, which addresses many of the spe cial difficulties 
posed by anomalous events ( ^ackett 19991:| ). 

9. CONCLUSION 

In this paper, we have presented and tested a method to 
calculate the detection efficiency of microlensing datasets 
to stellar and planetary lensing companions. This method 
is conceptually simple, direct, and not excessively time 
consuming. The final result, e{b,q), conveniently summa- 
rizes the efficiency of a given dataset in revealing binary 
lensing systems with any mass ratio q and dimensionless 
separation b. This efficiency can then be used to evaluate 
the likelihood that any given model of planetary systems 
would give rise to the observed dataset. 

The algorithm that we have presented has several ad- 
vantages. First, by working directly with measured light 
curves, it makes few assumptions, and explicitly takes into 
account the dependence of the efficiency on the sampling, 
photometric precision, monitoring duration, and impact 
parameter of the events. Second, it is computationally 
efficient, involving direct integration over only one un- 
known parameter, the angle of the trajectory. Finally, it 
is convenient, incorporating all information into one two- 
dimensional function, e{b,q), which can then be used to 
evaluate the likelihood of any model of planetary systems. 

Some caveats must be noted. First, the method assumes 
that the primary lensing parameters fo, Wmin, -Fo and 
fs are well-constrained by the light curve, while empha- 
sizing that e is poorly known in the absence of such con- 
straints. Possible altered sampling and photometric pre- 
cision of anomalous events may necessitate a resampling 
of the corresponding light curve or the assumption that 
the efficiency of anomalous events is equal to the aver- 
age efficiency of the normal events. Resampling is not an 
efficient use of all the data, and for individual events the 
average efficiency of non-anomalous events is likely to be a 
rather poor proxy for Cj. Finally, the algorithm implicitly 
assumes that the efficiency is independent of the plane- 
tary system model (e.g., the effect of multiple planets is 
ignored) . 

We have also applied our method to artificial light curves 
in order to explore the dependence of binary detection effi- 
ciency on the impact parameter of the event, the choice of 
the detection criterion, neglecting to fit comparison binary 
curves separately, finite source effects, and blending. Our 
conclusions are as follows: 

(1) Most of the efficiency of microlensing datasets to de- 

tecting planetary anomalies results from non-caustic 
crossing events. 

(2) Detection efficiency is strongly dependent on the im- 

pact parameter Mmin of the event. This implies 
that the integrated detection efhciency will depend 
strongly on the actual distribution of impact param- 
eters of the monitored events. 

(3) For mass ratios q > 10~^, the inferred detection effi- 

ciency e is robust to changes in the detection crite- 
rion Axthresh- smaller mass ratios q 



10 near 
can 



the detection limit, different choices of Axtin-ogh 
lead to differences in e of a factor of two. 
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(4) In order to obtain accurate detection efficiencies, the 

binary lens light curve must be fit separately unless 
the mass ratio is small {q~ 10""^) and the deviations 
are well above the detection limit (which may be dif- 
ficult to satisfy simultaneously). 

(5) Finite source effects are negligible for q~10~^, can 

either increase or decrease e for q ~ 10"'^, and can 
be devastating for qS 10~^, at least for the large di- 
mensionless source sizes {p~ 0.03) that are typical of 
bulge giants. 

(6) The detection efficiency is very sensitive to the frac- 

tion of blended light, primarily due to item (2) above; 
higher blend fractions imply smaller impact param- 
eters for detection, and smaller impact parameters 
have higher intrinsic detection efficiencies. This fact, 
combined with the degeneracy in the blend frac- 
tion that exists between poorly/inaccurately sam- 
pled blended light curves means that monitoring 
teams should make every effort to quantify the blend 



fractions in every monitored light curve. 

The final result of applying this method to an observed 
data set is the integrated detection efficiency, e{b, q). This 
result should be regarded as a primary outcome of mi- 
crolensing monitoring teams conducting planet searches; 
it incorporates all the information from the observations, 
and is essential for establishing any conclusions about the 
abundance and nature of Galactic planetary systems. 
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